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Abstract 

Background: Citrus variegated chlorosis (CVC), caused by Xylella fastidiosa, is one the most important citrus 
diseases, and affects all varieties of sweet orange [Citrus sinensis L Osb). On the other hand, among the Citrus 
genus there are different sources of resistance against X. fastidiosa. For these species identifying these defense 
genes could be an important step towards obtaining sweet orange resistant varieties through breeding or 
genetic engineering. To assess these genes we made use of mandarin (C. reticulata Blanco) that is known to be 
resistant to CVC and shares agronomical characteristics with sweet orange. Thus, we investigated the gene 
expression in Ponkan mandarin at one day after infection with X. fastidiosa, using RNA-seq. A set of genes 
considered key elements in the resistance was used to confirm its regulation in mandarin compared with the 
susceptible sweet orange. 

Results: Gene expression analysis of mock inoculated and infected tissues of Ponkan mandarin identified 667 
transcripts repressed and 724 significantly induced in the later. Among the induced transcripts, we identified genes 
encoding proteins similar to Pattern Recognition Receptors. Furthermore, many genes involved in secondary 
metabolism, biosynthesis and cell wall modification were upregulated as well as in synthesis of abscisic acid, 
jasmonic acid and auxin. 

Conclusions: This work demonstrated that the defense response to the perception of bacteria involves cell wall 
modification and activation of hormone pathways, which probably lead to the induction of other defense-related 
genes. We also hypothesized the induction of auxin-related genes indicates that resistant plants initially recognize 
X. fastidiosa as a necrotrophic pathogen. 
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Background 

The Brazilian citrus industry accounts for 30% of sweet 
orange production and 85% of exports of frozen-con- 
centrated orange juice in the world, despite the large 
number of pests and diseases that affect the Brazilian 
orchards. Among these diseases, Citrus Variegated Chlo- 
rosis (CVC), caused by the bacterium Xylella fastidiosa, 
costs around 120 million US dollars a year to chemically 
control the bacterial vectors and for replanting new or- 
chards [1]. 
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The symptoms of this disease are associated with the 
blockage of xylem vessels by X. fastidiosa biofilm, lea- 
ding to increased water stress and decreased nutrients in 
the diseased plant [2-4]. 

Citrus species show varying responses to CVC. While 
the sweet orange {Citrus sinensis L. Osb) is very suscep- 
tible, the Ponkan mandarin (Citrus reticulata Blanco) is 
considered resistant because it shows no symptoms, yet 
the bacteria can be isolated from the plants at 30 days 
after inoculation. However, after 60 days of inoculation 
the bacteria cannot be isolated from the plant. The re- 
sistance of mandarin is not related to the number and/ 
or diameter of xylem vessels, suggesting that resistance 
is caused by active defense responses [5]. Based on this, 
the pattern of gene expression in Ponkan mandarin was 
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assessed by sequencing expressed sequence tags in man- 
darins inoculated with X fastidiosa at 30 and 60 days 
after infection. The results revealed differential expres- 
sion patterns for several defense-related genes of the 
salicylic acid (SA), jasmonate (JA), and ethylene (ET) 
signaling pathways [6,7]. These results indicate a cross- 
talk between regulatory pathways that control different 
cellular processes in the mandarin-X. fastidiosa inter- 
action. However, it is unclear whether these pathways 
are activated during the initial response of Ponkan man- 
darin to this phytopathogen. Thus, the present study 
aimed to evaluate which genes are activated in the pre- 
liminary stages of infection, as this phase may involve an 
important strategy for avoiding pathogen establishment 
and colonization, and consequently the progress of the 
disease. Identifying these defense genes could be an im- 
portant step towards obtaining sweet orange resistant 
varieties through breeding or genetic engineering. 

Results and discussion 

Overview of RNA-seq analysis 

In recent years the number of works using global ex- 
pression analysis to study plant-pathogen interactions 
has grown considerably. By comparing specific mRNAs 
present in different tissues, such as infected or not 
infected, differentially expressed genes can be identified 
and their functions inferred. 

In the present study, we used RNA-seq to analyze the 
differential expression of Ponkan mandarin mRNAs one 
day after X. fastidiosa infection (compared with mock 
inoculated plants). The presence or absence of bacteria 
in the plants used in this analysis was confirmed by real- 
time quantitative PCR (RT-qPCR) (Additional file 1). 
Three biological replicates for each condition were se- 
lected for performing transcriptome analyses. 

RNA-seq generated 35,344,265 and 37,326,339 single 
end reads of 101 bp for the non-infected and infected li- 
braries, respectively. Tophat was used to align the reads 
to the Citrus Clementina reference genome, with ap- 
proximately 74% of success. The expression level of 
27,889 coding sequences was quantified using the soft- 
ware Cufflink. A total of 1,391 transcripts showed sig- 
nificant variation in expression by Cuffdiff analysis: 724 
were induced and 667 were repressed (P < 0.001) in the 
infected Ponkan mandarin (Additional file 2). 

The differentially expressed genes were categorized 
using Gene Ontology (GO). Based on similarity, trans- 
cripts were distributed into different categories of bio- 
logical processes (level 2). The prominent functional 
categories for both induced and repressed genes were re- 
lated to metabolic process, followed by cellular process, 
response stimulus, biological regulation, and localization 
(Figure 1A and B). The majority of induced genes are in- 
volved in the formation of secondary xylem, cell wall, 



lignin, hormone synthesis, and detoxification. On the 
other hand, the repressed genes are related to growth, 
cell wall degradation, cell wall loosening, cell differenti- 
ation, and development. A complete list of X. fastidiosa- 
modulated pathways for Ponkan mandarin data sets is 
shown in Additional file 2. 



Induction of a PAMP - Triggered Immunity (PTI) - like 
response in Ponkan mandarin 

In plants, Pattern Recognition Receptors (PRRs) perform 
the first level of microbial recognition, perceiving mo- 
lecular signatures characteristic of a whole class of 
microbes, which are termed Pathogen-Associated (or 
Microbe-Associated) Molecular Patterns (PAMPs or 
MAMPs) [8]. PAMP recognition leads to a chain of sig- 
naling events broadly referred to as general defense re- 
sponses in plants [9]. Stimulation of PRRs is a key step 
in the early stages of PTI. We identified two upregulated 
genes (Ciclevl0004108m and Ciclevl0014130m), which 
are similar to a leucine-rich repeat receptor-like protein 
kinase (LRR-RLK) in Arabidopsis (Table 1), belonging to 
the LRR-XII subfamily, of which EF-Tu and Flagellin 
sensing 2 receptors are members [10]. Upregulation of 
these genes suggests that PAMPs from X fastidiosa 
might be perceived by PRRs in Ponkan mandarin and 
could trigger the PTI response. After this recognition, 
the plant modulates the expression of defense genes as 
well as initiating metabolic rearrangements, and ulti- 
mately activates basal resistance to potential pathogens 
[11]. Another upregulated gene in Ponkan mandarin 
encodes a leucine-rich repeat receptor-like protein 
(RLP12), which is associated with different functions 
such as cell differentiation, plant growth, development, 
and mainly disease resistance [12-14] (Table 1). 

A gene encoding a coiled-coil motif nucleotide- 
binding site-leucine-rich repeat (CC-NBS-LRR) protein 
was upregulated in Ponkan mandarin challenged with X 
fastidiosa (Table 1). CC-NBS-LRR is a cytoplasmic re- 
ceptor normally involved in responses triggered after 
recognition of Avr proteins secreted by the pathogen. 
However, this bacterium does not have the type III se- 
cretion apparatus or Avr proteins [15]. This leads us to 
believe that this receptor may be involved in perception 
of damage-associated molecular patterns (DAMPs), be- 
cause the recognition of cytoplasmic danger signals de- 
pends on cytoplasmic sensors like NB-LRR resistance 
proteins [16-18]. This hypothesis is consistent with X. 
fastidiosa^ ability to produce danger molecules by de- 
grading plant cell walls [4]. However, this hypothesis 
needs to be further investigated. 

A main feature of the PTI response is strengthening of 
the cell wall [19]. Many genes related with cell wall 
modification were upregulated in Ponkan mandarin 
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Figure 1 Categorization of genes induced or repressed in Ponkan mandarin infected with X. fastidiosa. The genes were 
automatically classified based on GO (Gene Ontology) classifications and are grouped according to their functions in plant cells. A. Induced genes. 
B. Repressed genes. 



infected with X fastidiosa, which reinforced our hypo- 
thesis that PTI is involved in this early response. 

Genes related to secondary metabolism and the cell wall 

In this study, we observed a significant change in expres- 
sion of genes involved in secondary metabolism, and cell 
wall biosynthesis and modification in Ponkan mandarin 
infected with X fastidiosa. These genes were mapped 
using MapMan to generate a representative overview 
(Figure 2 and Additional file 3). 

In response to X fastidiosa infection, Ponkan manda- 
rin appears to activate genes for the formation and alter- 
ation of secondary xylem cells as a defense mechanism. 
Induction of the transcription factor AS1, which has a 
MYB domain (Ciclevl0012089m) was also observed. It 
is suggested that the abundance of MYB proteins in the 
xylem could be involved in transcriptional regulation of 



the formation of the secondary xylem [20] (Table 1). 
Furthermore, a myo gene (Ciclevl0010780m), which 
encodes actin, was strongly induced (Table 1). Many 
studies suggest that actin displays a similar expression 
pattern to microtubules, which are suggested to deter- 
mine the location for deposition of the secondary cell 
wall, cellulose, lignin, hemicellulose, and proteins [21]. 
Another transcription factor induced in this study was 
MYB66 (Ciclevl0017556m), which is probably involved 
in the regulation of flavonoids biosynthesis and lignifi- 
cation [22-24] (Table 1). To corroborate these observa- 
tions, we verified the activation of genes related to 
phenylpropanoid and flavonoids biosynthesis, which act 
as anti-pathogenic molecules [25]. Some of them such as 
phenylalanine ammonium lyase are involved in lignin 
biosynthesis (Table 1; Figure 2 and Additional file 3). In 
addition to lignification, callose deposition is also an 
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Table 1 Differentially expressed genes in Ponkan mandarin in response to infection by X. fastidiosa 



Gene symbol 


gene \d_Citrus 
Clementina* 


AGI** 


Fold change*** 


P< 0.001 


Gene description 


LRR-RLK 


Ciclevl 00041 08m 


AT4G08850 


1.05611 


0.000755101 


leucine-rich repeat family protein 


LRR-RLK 


Ciclevl 00141 30m 


AT3G47570 


2.64522 


3.84525e-09 


leucine-rich repeat protein kinase family protein 


RLP12 


Ciclevl 0003540m 


AT1G71400 


1.20144 


2.69007e-06 


receptor like protein 12 


CC-NBS-LRR 


Ciclevl 0007304m 


AT4G27190 


1.29361 


4.28293e-06 


nb-arc domain-containing disease resistance protein 


MYB domain 


Ciclevl 001 2089m 


AT2G37630 


1 .35348 


0.000930383 


asl (asymmetric leaves 1); transcription factor 


MYO 


Ciclevl 001 0780m 


AT4G33200 


2.76714 


0.000191443 


myosin, putative 


MYB66 


Ciclevl 001 7556m 


AT5G 14750 


1 .05408 


1 .04467e-07 


myb domain protein 66); transcription factor 


PAL 


Ciclevl 0030821m 


AT2G37040 


0.91583 


2.53552e-07 


phenylalanine ammonia-lyase 


GSL07 


Ciclevl 0030560m 


AT1G06490 


1.0981 


2.36835e-05 


glucan synthase-like 7 


CSLA09 


Ciclevl 0031 284m 


AT5G03760 


1.33401 


1 .33936e-07 


cellulose synthase like A9 


CESA8 


Ciclevl 00141 55m 


AT4G 18780 


1 .27949 


6.74384e-08 


cellulose synthase 


CESA4 


Ciclevl 001 8639m 


AT5G44030 


1.15575 


7.36883e-05 


cellulose synthase 


ATFXG1 


Ciclevl 0031 741m 


AT1G67830 


-0.940316 


1 .02739e-05 


alpha-L-fucosidase/ carboxylesterase 


# 


Ciclevl 0020423m 


AT4G24780 


-1.03748 


3.74844e-06 


pectate lyase family protein 


# 


Ciclevl 0030947m 


AT1G70370 


-1.27505 


1 .32952e-06 


BURP domain-containing protei 


QRT2 


Ciclevl 000531 7m 


AT3G07970 


-1.57918 


4.6342 1e-06 


QRT2; polygalacturonase 


# 


Ciclevl 0004783m 


AT1G04680 


-1.68649 


1 .5906e-05 


pectate lyase family protein 


LOX 


Ciclevl 0014202m 


AT3G45140 


0.99318 


0.000466003 


lipoxygenase 2 


NCED6 


Ciclevl 0006710m 


AT3G24220 


2.69285 


1 .90022e-09 


nine-cis-epoxycarotenoid dioxygenase 6 


CCD7 


Ciclevl 0027500m 


AT2G44990 


0.93434 


0.00201556 


carotenoid clevage dioxygenase 7 


AP2 


Ciclevl 001 0403m 


AT1G51 120 


2.21376 


5.89723e-07 


ap2 domain-containing transcription factor, putative 


AIP2 


Ciclevl 0008969m 


AT5G20910 


-1.33947 


4.99589e-07 


zinc finger (c3hc4-type ring finger) family protein 


IAA9 


Ciclevl 0025873m 


AT5G65670 


2.26201 


0.0033 


indole-3-acetic acid inducible 9 


E3 RING 


Ciclevl 0005027m 


AT2G30580 


1 .46874 


1 .05377e-05 


ubiquitin-protein ligase 


HECT (UPL5) 


Ciclevl 001421 3m 


AT4G 12570 


1.41737 


6.03935e-06 


ubiquitin protein ligase 5 


E3/SCF/FB0X (LKP1) 


Ciclevl 0007762m 


AT5G57360 


0.98971 


0.00243897 


ubiquitin-protein ligase 


ARF-GAP DOMAIN 8 


Ciclevl 001 1983m 


AT4G 17890 


0.87335 


0.00250141 


arf gtpase activator 


ARF8 


Ciclevl 0014200m 


AT5G37020 


0.83349 


0.00155566 


transcription factor 


ARF19 


Ciclevl 0007286m 


AT1G 19220 


1.17169 


0.000240451 


transcription factor 


BIG 


Ciclevl 001 0885m 


AT3G02260 


0.62707 


0.000252397 


auxin transport protein (big) 


XTH16 


Ciclevl 0005561m 


AT3G23730 


-0,983574 


5.04445e-07 


xyloglucan endotransglycosylase, putative 


EXGT-A4 


Ciclevl 002885 1m 


AT5G 13870 


-1,37398 


5.49438e-1 1 


engoxyloglucan transferase A4 


ATEXPA4 


Ciclev10012518m 


AT2G39700 


-1,99214 


0.000463125 


expansin A4 



* Citrus Clementina transcripts identification number (reference genome used for mapping the reads) - http://www.phytozome.org/search.php. 

** Identification number of the Arabidopsis thaliana ortholog of up and down-regulated citrus gene in response to X. fastidiosa infection (The Arabidopsis 
Genome Initiative). 

*** Log 2 fold change values (P< 0.001) obtained from of each infected sample compared to mock-inoculated control. 

# unnamed gene symbol. 

Up and downregulated genes in Ponkan mandarin 24 h after infection with X. fastidiosa compared with the control. 



important defense mechanism in plants and one callose 
synthase (Ciclevl0030560m) was induced in Ponkan 
mandarin infected with X. fastidiosa (Table 1; Figure 2 
and Additional file 3). Additionally, genes encoding cel- 
lulose synthases (Ciclevl003 1284m, Ciclevl0014155m 
and Ciclevl 00 18639m) were significantly induced, such 
as CESA8 and CESA4, which are key enzymes in the 



biosynthesis of the xylem cell wall [26]. Among the re- 
pressed genes it is remarkable the presence of those en- 
coding proteins related to cell wall degradation (Table 1; 
Figure 2 and Additional file 3). These results indicate 
that the molecular defense response of Ponkan mandarin 
against X. fastidiosa involves the participation of genes 
related to cell wall biosynthesis. This could represent an 
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Figure 2 Responses related to secondary metabolism and the cell wall in Ponkan mandarin 1 day after infection with X. fastidiosa. 

Log 2 fold change of gene expression {X. fastidiosa versus mock inoculated control) was analyzed by MapMan. Blue squares represent upregulated 
genes, red squares represent downregulated genes and grey circles represent not differentially regulated genes. The color saturation indicates 
fold change >0.8 and < -0.8. GSL07 (Ciclevl 0030560m), glucan synthase-like 7; CSLA09 (Ciclevl 0031 284m), cellulose synthase like A9; CESA8 
(Ciclev10014155m), cellulose synthase; CESA4 (Ciclevl 001 8639m), cellulose synthase. 



important strategy of the plant for restrict the movement 
of X. fastidiosa through the xylem cells. 

Hormone related-genes in the Ponkan mandarin defense 
response 

After pathogen recognition, plants transmit signals to 
activate defense responses. This transmission can be 
performed by secondary messengers, such as G-proteins, 
Ca 2+ , reactive oxygen species (ROS), nitric oxide, and 
hormones [27,28]. Genes associated with calcium 
signaling (Ciclevl0000095m, Ciclevl0014823m and 
Ciclevl0032932m), G-proteins (Ciclevl0014139m, 
Ciclevl0007246m and Ciclevl0020496m) and ROS 
(Ciclevl0026073m and Ciclevl0026072m) were all 
upregulated in infected plants (Table 1 and Table 2). 
In addition, hormone related-genes induced in Ponkan 
mandarin infected with X. fastidiosa were associated with 
JA, ABA and auxin. JA and ABA pathways are important 
regulators of expression of defense genes and have been 
identified downstream in PTI [29,30]. The gene encoding 
lipoxygenase (LOX), a key enzyme in the synthesis of JA 
[31], was induced in Ponkan mandarin infected by X. 
fastidiosa (Table 1; Figure 3 and Additional file 4). In 
addition to the defense response, this hormone activates 
secondary metabolism in the plant in response to a variety 
of biotic and abiotic stresses [31,32]. Interestingly, LOX 
was also upregulated in Ponkan mandarin in later stage of 
X fastidiosa infection [6,7]. These observations highlight 
the importance of the JA pathway during the defense 
response. 

With regard to ABA-related genes, AP2 (CiclevlOO 
10403m), nced6 (Ciclevl0006710m), and ccaV (Ciclev 
10027500m) were induced in response to X fastidiosa 
infection (Table 1; Figure 3 and Additional file 4). AP2 is in- 
volved in the activation of genes related to ABA 



biosynthesis [33], while nced6 and ccdl are associated 
with biosynthesis and transport of this hormone, respect- 
ively. Additionally the aipl gene (Ciclevl0008969m) was 
repressed (Table 1; Figure 3 and Additional file 4) and 
this gene encodes an E3 ligase that negatively regulates 
ABA signaling by targeting ABB, a central regulator of 
this pathway, for degradation [34] . 

So the repression of this gene strengthens the idea that 
ABA biosynthesis is activated in Ponkan mandarin in re- 
sponse to X. fastidiosa infection. 

Other hormone related-genes upregulated in Ponkan 
mandarin in response to infection by X. fastidiosa were 
associated with the auxin signaling pathway (Table 1; 
Figure 4 and Additional file 5). Indole-3-acetic acid 
(IAA) is the main auxin in plants, controlling many im- 
portant physiological processes, including cell growth 
and division, tissue differentiation and response to light 
[35,36]. In addition, auxin is also associated with in- 
creased susceptibility to biotrophic microorganisms, be- 
cause it promotes loosening of the cell wall and thus 
potentiates pathogen growth [30,37]. Many bacteria pro- 
duce IAA as a strategy to interfere with the plant auxin 
pathway to facilitate their infection [38]. However, our 
results showed the induction of several genes involved in 
the activation of the auxin signaling pathway in a resis- 
tant plant after infection, suggesting that the plants do 
not recognize X. fastidiosa as a biotrophic pathogen. 
The map locations of modulated auxin genes found in 
this study are shown in representative schematics of 
auxin synthesis and degradation (Figure 4). At high 
concentrations, auxin promotes an association between 
Auxin/IAA (Aux/IAA) and an F-box protein known as 
transport inhibitor response 1 (TIR1). This complex ac- 
tivates the E3 ligase that leads to degradation of Aux/ 
IAA, allowing the release of auxin responsive factors 
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Table 2 Defense-genes upregulated in Ponkan mandarin in response to X. fastidiosa infection 



Gene symbol Gene \d_Citrus AGP* Fold change*** P < 0.001 Gene description 

Clementina* 



# 


Ciclevl 0000095m 


AT1G20760 


0.872938 


0.000738646 


calcium-binding EF hand family protein 


CPK21 


Ciclevl 0014823m 


AT4G04720 


0.731932 


0.0054887 


calmodulin-dependent protein kinase 


# 


Ciclevl 0032932m 


AT1G64850 


0.600133 


0.0256525 


calcium-binding EF hand family protein 


# 


Ciclevl 00141 39m 


AT5G46070 


0.882012 


0.046128 


GTP binding / GTPase 


ATMIN7 


Ciclevl 0007246m 


AT3G43300 


0.625071 


0.030625 


guanyl-nucleotide exchange factor/ protein binding 


# 


Ciclevl 0020496m 


AT5G50120 


0.617877 


0.00757447 


WD-40 repeat family protein 


PER17 


Ciclevl 0026073m 


AT2G22420 


0.88525 


0.0376681 


peroxidase 17 


# 


Ciclevl 0026072m 


AT4G33420 


0.559908 


0.00032072 


peroxidase, putative 


MAP9 


Ciclevl 0021 170m 


AT1G73500 


0.567652 


0.000807898 


MAP kinase kinase 9 


EMB3004 


Ciclevl 0000874m 


AT3G06350 


1 .0645 


0.00113401 


shikimate dehydrogenase, putative 


HSP70 


Ciclevl 0027981m 


AT3G 12580 


1 .2447 


4.22995e-12 


heat shock protein 70 


HSP17 


Ciclevl 0009756m 


AT5G 12020 


1 .62260 


0.00 


17.6 kda class II heat shock protein 


HSP90 


Ciclevl 0004456m 


AT5G52640 


1.63013 


9.87432e-13 


heat shock protein 90.1 


DNAJ 


Ciclevl 0028792m 


AT4G39150 


1.01887 


0.00183666 


dnaJ heat shock N-terminal domain-containing protein 


HPT1 


Ciclevl 0008336m 


AT2G 18950 


2.40936 


2.24863e-10 


homogentisate phytyltransferase 1 


P450 


Ciclevl 0014861m 


AT4G 12320 


1 .56043 


0.000237455 


cytochrome P450, family 706, subfamily A, polypeptide 6 


TAU8 


Ciclevl 0005809m 


AT3G09270 


1 .86905 


1. 54321 e-08 


glutathione S-transferase TAU 8 


UDP-glycosyltransferase 


Ciclevl 0025462m 


AT3G50740 


1.139 


1 .72708e-07 


udp-glucosyl transferase 72E1 


CHIA 


Ciclevl 0005491m 


AT5G24090 


1 .02503 


0.00135025 


chitinase A 


TRYPSIN 


Ciclevl 0022001m 


AT1G 17860 


1.32418 


2.00905e-06 


trypsin and protease inhibitor family protein 



* Citrus Clementina transcripts identification number (reference genome used for mapping the reads) - http://www.phytozome.org/search.php 

** Identification number of the Arabidopsis thaliana ortholog of up and down-regulated citrus gene in response to X. fastidiosa infection (The Arabidopsis 
Genome Initiative). 

*** Log 2 fold change values (P< 0.001) obtained from of each infected sample compared to mock-inoculated control. 

# unnamed gene symbol. 

Upregulated genes in Ponkan mandarin 24 h after infection with X. fastidiosa compared with the control. 



(ARFs) from the complex [39]. A gene encoding Aux/IAA 
(IAA9 - Ciclevl0025873m) was the most upregulated in 
our analysis. RING E3 subunit (Ciclevl0005027m), a 
HECT subunit (UPL5 - Ciclevl0014213m) and E3/SCF/ 
FBOX (LKP1-Ciclevl0007762m) genes, which are all part 
of the proteasome complex, were also induced. Addition- 
ally, an ubiquitin-specific protease (ARF-GAP DOMAIN 
8 - CiclevlOOl 1983m) was also significantly induced in 
the infected plant (Table 1; Figure 4 and Additional file 5). 
This enzyme is highly conserved in eukaryotes, and plays 
a critical role by cleaving ubiquitinated proteins [40]. Fur- 
thermore, two ARFs (ARF19 and ARF8) related-genes and 
a gene encoding BIG, which is involved in polar auxin 
transport and has an essential function in auxin signaling 
[41] were upregulated (Table 1). 

Many studies have reported that auxin promotes sus- 
ceptibility to bacterial diseases [42,43]. However, this af- 
firmation is only true for biotrophic organisms: auxin 
signaling is an important component involved in plant 
resistance to necrotrophic pathogens [30,44]. Auxin sig- 
naling is upregulated in Ponkan mandarin in response to 
X. fastidiosa and this plant is resistant to this pathogen; 



therefore, we speculated that in the early stage of infec- 
tion Ponkan mandarin recognizes X fastidiosa as a 
necrotrophic organism, even though this bacterium does 
not cause massive destruction of host tissue. This is con- 
sistent with the direct injection of the bacteria by the in- 
sect vector into the xylem vessels, which is composed 
mainly of dead cells, and the fact that X fastidiosa is 
able to degrade plant cell wall, which is an important 
factor for its colonization in susceptible plants [4]. 
Therefore, substantial tissue destruction is not necessary 
to trigger necrotroph-related responses mediated by X 
fastidiosa in the resistant plant host. 

The evidence in this paper indicates that activation of 
the auxin signaling pathway does not promoting suscep- 
tibility of Ponkan mandarin after infection with the bac- 
terium. The involvement of this hormone in pathogen 
susceptibility disease development appears to include 
rapid elongation of plant tissues by increasing the exten- 
sibility of the cell wall [45]. Proteins that participate in 
acid-induced cell wall extension are endo-(3-l,4-glucanases 
(EGases), xyloglucan endotransglycosylases (XETs), and 
expansins [46]. We did not observed induction of a major 
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Figure 3 Stress responses in Ponkan mandarin 1 day after infection with X. fastidiosa. Log 2 fold change of gene expression (X fostidiosa 
versus mock- inoculated control) was analyzed by MapMan. Blue squares represent upregulated genes, red squares represent downregulated 
genes and grey circles represent not differentially regulated genes. The color saturation indicates fold change >0.8 and < -0.8. 



set of these proteins. In fact some of them such as XETs 
(Ciclevl0005561m and Ciclevl0028851m) and expansin 
(Ciclevl0012518m) were repressed in Ponkan mandarin 
in response to X. fastidiosa, suggesting that auxin signal- 
ing related genes found in this work are not primarily pro- 
moting cell expansion (Table 1). On the other hand, auxin 
synergistically with jasmonic acid/ethylene-signaling is 
required for necrotrophic resistance. In addition some 
members of ARF positively regulate camalexin biosyn- 
thesis resulting in resistance to necrotrophic pathogens 
[47]. This defense system, involving the interaction be- 
tween auxin and JA signaling pathways, which probably 
occurs in Ponkan mandarin, given the significant induc- 
tion of genes associated with both signaling pathways. 
Camalexin is produced through the tryptophan pathway 
and in our analysis two genes involved in the biosynthesis 
of tryptophan, dehydroquinate-shikimate dehydrogenase 
(Ciclevl0000874m) and MAPK 9 (Ciclevl0021170m), 
were induced in Ponkan mandarin, suggesting that ca- 
malexin biosynthesis may be induced in response to X. 
fastidiosa infection. Taken together all these evidences 
suggest us that the resistant Ponkan mandarin recognizes 
X, fastidiosa in early stage of infection as a necrotrophic 
pathogen. 

To confirm that auxin signaling related genes are in- 
deed upregulated only in Ponkan mandarin in response 
to X. fastidiosa infection we also evaluated the expres- 
sion of auxin marker-genes by RT-qPCR in Pera sweet 
orange susceptible variety and in Ponkan mandarin at one 
day after X fastidiosa infection. As shown in Figure 5, all 
auxin related-genes were significantly induced only in 



Ponkan mandarin. In Pera sweet orange, the genes were 
significantly repressed (E3 and ARF 19) or showed no sig- 
nificant change (IAA9, TIR1 and BIG). This result evi- 
dences that auxin is induced as a resistance response 
against X. fastidiosa during the early stage of infection. 
Recognition of PAMPs or DAMPs that somehow resem- 
ble necrotrophic pathogens may mediate this response. 
However, this recognition occurs mainly during the early 
stage of infection since we observed a gradual decrease in 
expression of auxin related-genes along the time course of 
infection (Figure 6). After 21 days, no auxin related-gene 
was expressed, whereas expression of salicylic acid (SA) 
marker-gene increased (Figure 6). This result agrees with 
De Souza et al. [6,7] where an upregulation of SA related- 
genes was observed in Ponkan mandarin at 30 days after 
X. fastidiosa inoculation. After this time point, the bacter- 
ial population decreases to a point where it could not be 
isolated [3]. These results suggest that the resistant plant 
changes its mechanism of defense during X. fastidiosa in- 
fection: the initial response involves the participation of 
auxin while later on SA becomes important. It is to note 
that the change occurs approximately at the time when X 
fastidiosa forms a structured biofilm. In this growth con- 
dition this bacterium expresses specific genes and proteins 
necessary for its adaptation and pathogenicity in the host 
[3,4,48,49]. Therefore other proteins expressed in biofilm 
condition could be later recognized by the plant. Ne- 
vertheless, how the resistant plant indeed recognizes X 
fastidiosa to trigger different pathways in the resistance 
response remains to be discovered. Other downstream 
defense-genes upregulated in Ponkan mandarin after X 
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Figure 4 The ubiquitin dependent degradation system in Ponkan mandarin 1 day after infection with X. fastidiosa. Activation of the 
auxin pathway in Ponkan mandarin after infection with the bacteria. Signaling of this hormone happens when, at high auxin concentrations, the 
transcription factor Aux /IAA9 is sent to ubiquitin dependent degradation, which leads to degradation of this factor, allowing the release of auxin 
responsive factors (ARFs) from the complex. Log 2 fold change of gene expression (X. fastidiosa versus mock-inoculated control) was analyzed by 
MapMan. Blue squares represent upregulated genes, red squares represent downregulated genes and grey circles represent not differentially 
regulated genes. The color saturation indicates fold change >0.8 and < -0.8. Aux /IAA9, transcription factor Auxin/lndole-3 acetic acid 9; DUB, 
Deubiquitinating enzyme; El, Ubiquitin-activating enzymes; UBA1, ubiquitin-activating enzyme E1; E2, Ubiquitin-conjugating enzymes; UBC2, 
ubiquitin-conjugating enzyme E2; RING, C3HC4 RING-domain-containing ubiquitin E3 ligase; HECT, HECTtype E3; APC, anaphase-promoting 
complex; E3, Ubiquitin ligases; Cullin, SKP, FBOX and RBX, subunit of the E3 ligase; ARF, auxin responsive factors. 



fastidiosa infection are represented at the Table 2. These 
genes might contribute to increase the resistance response 
in Ponkan mandarin to X. fastidiosa culminating in its 
elimination in the plant. 

To confirm the participation of different genes related 
to pathogen recognition, cell wall synthesis, and hor- 
mone signaling pathways in the mandarin resistance 
response, we also tested their expression in Pera sweet 
orange. The analysis confirmed that genes encoding 
LRR-RLK and CC-NBS-LRR (pathogen recognition), AP2 
(ABA signaling), MYO and CESA4 (cell wall synthesis) 



were not only specifically induced in mandarin but also 
repressed in sweet orange (Additional file 6). 

Validation of RNA-seq data by RT-qPCR 

RT-qPCR was used to validate the RNA-seq data. 
Twelve genes involved in different biological processes 
were selected (Additional file 7). Similar expression pat- 
terns were observed for all genes evaluated by both tech- 
niques (Additional file 8). Additionally, a high Spearman's 
rho value (0.88) indicated a good correlation between the 
fold change from RNA-seq experiments and RT-qPCR. 
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■ Ponkan mandarin 
Pera sweet orange 
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IAA9 TIR1 E3 ARF19 BIG 

Figure 5 Relative quantification of genes related auxin pathway in Ponkan mandarin and sweet orange by RT-qPCR. cDNA samples 
were prepared using RNA from xylem tissue from Ponkan mandarin and sweet orange, after 1 day of infection with X. fastidiosa or not (control) 
(three biological replicates). The bars indicate the standard deviation of the means. (*) indicates significant difference (P > 0.05) between the 

mean values obtained for each gene compared with the control. 
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IAA9 pr1 

Figure 6 Relative quantification of genes encoding IAA9 and PR1 in Citrus plants infected with X. fastidiosa by RT-qPCR. cDNA samples 
were prepared using RNA extracted from a mixture of leaves and petioles of Ponkan mandarin and Pera sweet orange at 1, 7, 14, and 21 days 
after inoculation (d.a.i.) of X. fastidiosa or mock inoculated (control). The experiment was conducted with three biological replicates. The bars 
indicate the standard deviations of the means. (*) indicates significant difference (P > 0.05) between the mean values obtained for each gene 
compared to control. 



These results confirmed the reliability and accuracy of the 
RNA-Seq data in this study (Figure 7). 

Conclusions 

This work demonstrated that the defense response of 
Ponkan mandarin to X. fastidiosa involves induction of 
genes related to PAMP receptors, cell wall synthesis, and 
ABA, JA and auxin signaling pathways, which will prob- 
ably culminate in the induction of other defense-related 
genes. We also hypothesize that Ponkan mandarin ini- 
tially recognizes X. fastidiosa as a necrotrophic pa- 
thogen, inducing genes related to the auxin signaling 
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Figure 7 Correlation between RNA-seq and RT-qPCR data. 

Twelve differentially expressed genes in Ponkan mandarin 
challenged with X. fastidiosa were selected for validation. 
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pathway and later it changes the response to a bio- 
trophic pathogen. 

The genes found in this work are important tools to 
be used in breeding programs to accelerate the selection 
of hybrid from C. reticulata and C. sinensis carrying 
defense responses from mandarin. Moreover these genes 
can be also used for genetic engineering aiming the 
development of resistant sweet orange varieties. 

Methods 

Plant materials and bacterial detection 

The experiments were conducted with physiologically 
mature plants, of uniform size and production of young 
leaves. When the shoots reached a length of approxi- 
mately 10 cm, they were artificially inoculated by needle 
prick with 10 uL suspension (10 8 cells ml/ 1 ) of X. fas- 
tidiosa strain 9a5c in PBS buffer at five different points 
on the same stem. Negative controls comprising Ponkan 
mandarin (C. reticulata Blanco) and Pera sweet orange 
{Citrus sinensis L. Osb) were mock inoculated with PBS 
buffer. 

Total genomic DNA (plant + bacteria) was extracted 
from inoculated cambial tissue enriched with xylem after 
1 day using a CTAB method adapted from [50]. These 
samples were used for detection of bacteria using real- 
time PCR reactions. The analyses were carried out with 
an ABI PRISM 7500 Sequence Detector System (Applied 
Biosystems, Foster City, CA, USA). The reaction was per- 
formed in a total volume of 25 uL, containing 12.5 uL of 
TaqMan PCR Master Mix fast (Applied Biosystems), 
200 ng of DNA template and 525 nM of primers CVC-1 
and CCSM-1 [51]. Each sample was tested in triplicate 
and with five biological replicates. Negative (no template 
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DNA) and positive (DNA from X. fastidiosa) controls 
were included in all experiments to exclude or detect any 
possible contamination. The samples were considered 
positive for the presence of X fastidiosa when they pre- 
sented Ct (Cycle Threshold) below or equal to 37; negative 
samples did not present amplification up to this Ct value 
(unpublished data). From these results, we selected the 
three biological replicates for different conditions and 
performed transcriptome analyses. 

RNA isolation and expression analysis in RNA-seq 

For transcriptome analysis we used cambial tissue en- 
riched with xylem from Ponkan mandarin. Total RNA 
was extracted with Trizol (Invitrogen Life Technologies 
Foster City, CA, USA) and treated with DNase RNase- 
Free Set (Qiagen, Valencia, CA, USA), according to the 
manufacturers instructions. We extracted RNA samples 
from a pool of three independent biological replicates 
infected after one day and their respective controls (not 
infected). The concentration of RNA was measured in a 
NanoDrop ND-1000 spectrophotometer (NanoDrop 
Technologies, Wilmington, DE). RNA quality was evalu- 
ated using an Agilent Bioanalyzer Model 2100 (Agilent 
Technologies, Palo Alto, CA). 

A total of 10 \ig of RNA from Ponkan mandarin 
(control and infected with X. fastidiosa) were sent to 
Macrogen Inc. (South Korea) for sequencing using the 
Genome Analyzer IIx platform (Illumina Inc.). All proce- 
dures were performed according to Aluminas protocols. 
Purified cDNA libraries were dispersed onto an Illumina 
single-end flow cell composed of eight lanes using the 
Illumina Cluster Station (Illumina Inc.). One lane was 
used per sample of the treated and control plants. The 
101 bp reads were collected using the Illumina GA II 
and sequencing-by-synthesis technology. The sequences 
of Ponkan mandarin were mapped against the Citrus 
Clementina reference genome (http://www.phytozome. 
net/clementine.php) using TopHat [52]. After alignment, 
the relative abundance of the transcripts was measured 
with the Cufflink software, which measures the tran- 
scripts abundance as RPKM (Reads Per Kilobase of exon 
model per Million mapped reads). The differential ex- 
pression between Ponkan mandarins inoculated or not 
with bacteria, and its significance, was calculated in 
Cuffdiff [53]. The differentially expressed transcripts 
were annotated and automatically categorized using GO 
(Gene Ontology - http://www.blast2go.com/b2ghome). 
These sequences were also used to search for similar 
protein sequences available in GenBank using the 
BLASTX tool. 

In addition, the differentially expressed genes were also 
functionally analyzed using the MapMan software, which 
is a user-driven tool that displays large genomics 



datasets onto diagrams of metabolic pathways or other 
processes [54]. 

Expression analysis by RT-qPCR 

Twelve genes that were identified by RNA-seq to be in- 
duced or repressed in Ponkan mandarin in response to 
X. fastidiosa, were selected for validation by RT-qPCR: 
ATEXPA4; CLV1; CC-NBS-LRR; RLK; P12; LOX; AIP; 
MYO; AP2; HSP90; CCR4 and IAA9 (Additional file 7). 

Furthermore, we used RT-qPCR to compare the level 
of expression of some genes involved in the auxin path- 
way, pathogen recognition, ABA signal transduction, 
and cell wall synthesis in Ponkan mandarin and in Pera 
sweet orange, a susceptible variety, one day after infec- 
tion with the bacteria. We evaluated the IAA9 (Aux/ 
IAA), ARF19, TIR1, BIG and E3 genes (auxin pathway); 
LRR-RLK and CC-NBS-LRR (pathogen recognition); 
AP2 (ABA); MYO and CESA4 (cell wall synthesis). 

Additionally, we checked the relative quantification of 
genes encoding IAA9 and PR1 in Citrus plants, by RT- 
qPCR, using RNA extracted from a mixture of leaves 
and petioles of Ponkan mandarin and Pera sweet orange 
at 1, 7, 14, and 21 days after inoculation of X. fastidiosa 
or mock inoculated (control). 

The primers for these genes were designed using 
PrimerExpress software (Applied Biosystems, Foster City, 
CA, USA) (Additional file 9). 

The specificity of the primers was checked in silico 
against the NCBI database (http://www.ncbi.nlm.nih.gov/) 
using the Primer-BLAST tool. All the primer sequences 
showed specificity with the sequences of target genes. 
Additionally, we checked the pattern of dissociation ob- 
tained after RT-qPCR, using a meting curve for each 
primer. This showed a single peak for all the evalua- 
ted genes, confirming the existence of only one amplicon 
(Additional file 10). 

The efficiency of the primers was estimated in each 
experiment using the software Miner (http://www.miner. 
ewindup.info/). This software quantifies the results of 
RT-qPCR based on the kinetics of the PCR amplification 
individually for each sample, without the need for a 
standard curve. This allows a direct calculation of the 
efficiency and values of cycle quantification (cq) [55]. 
All primers showed amplification efficiencies between 
90 - 100% (Additional file 10). 

To find a reference gene to normalize the RT-qPCR 
results, the stability of five endogenous control genes in 
Citrus was analyzed to confirm their stability using 
geNorm software [56] and to ensure the existence of 
gene expression variation due to the experimental condi- 
tions. The primers for these genes were obtained from a 
previous work [57]. In this evaluation we used samples 
of Pera sweet orange and Ponkan mandarin (control 
and infected with X. fastidiosa), Ubiquitin (UBQ) and 
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degradation probably involved in auxin pathway in Ponkan mandarin 
24 h after infection with X. fastidiosa. 

Additional file 6: Relative quantification of genes related to 
pathogen recognition, cell wall synthesis, and hormone signaling 
pathways in Citrus plants infected with X. fastidiosa by RT-qPCR. 

cDNA samples were prepared using RNA from xylem tissue from Ponkan 
mandarin and Pera sweet orange, after 1 day of infection with or without 
(control) X. fastidiosa (three biological replicates). The bars indicate the 
standard deviation of the means. (*) indicates significant difference 
(P > 0.05) between the mean values obtained for each gene [LRR-RLK 
and CC-NBS-LRR (pathogen recognition); AP2 (ABA); MYO and CESA4 
(cell wall synthesis)] compared with the control. 

Additional file 7: Genes involved in different biological processes 
selected from RNA-seq Ponkan mandarin infected with X. fastidiosa. 

Twelve differentially expressed genes in Ponkan mandarin infected 
with X. fastidiosa were selected to validate the results obtained from 
RNA-seq analysis. 

Additional file 8: Validation of 12 differentially expressed genes 
selected from RNA-seq analysis by RT-qPCR. The fold changes are 
shown for 12 differentially expressed genes identified using RNA-seq 
compared to those obtained by RT-qPCR. For this, cDNAs were prepared 
from RNA of Ponkan mandarin xylem tissue infected with X. fastidiosa or 
not (control) after one day, with three biological replicates. RT-qPCR data 
were normalized to the two most stable endogenous control genes 
(UBQ and CYP). 

Additional file 9: Genes and primers used in RT-qPCR. Includes 
sequences of oligonucleotide primers used for RT-qPCR analysis. 

Additional file 10: Validation of the specificity and amplification 
efficiency of the RT-qPCR primers. Amplification of cDNA for Ponkan 
mandarin genes ATEXPA4, C/V1, CC-NBS-LRR, LRR-RLK, PI 2, LOX, AIP, MYO, 
API, H5P90, CCR4, IAA9, ARF19, TIR], BIG, E3, PR] and CESA4. (A) 
Verification the dissociation pattern obtained after RT-qPCR. Each gene 
showed a single peak after melting curve analysis, confirming the 
specificity the primers. (B) Efficiency of amplification using primers 
obtained through Miner software after RT-qPCR. 

Additional file 11: Average expression stability values (M) of the 
five citrus endogenous control genes calculated by geNorm. 

Expression stability values were calculated for samples from citrus 
genotypes infected with X. fastidiosa and controls (mock). A lower M 
value indicates more stable expression. 



cyclophilin (CYC) were the most stable and were se- 
lected for further analysis. However, the other three 
genes, eukaryotic translation elongation factor 2, NADP- 
isocitrate dehydrogenase and tubulin also showed sa- 
tisfactory mean values (M-value) (Additional file 11). 
These M-values are within acceptable values at a cutoff 
value of 0.15 [57]. 

For the analyses of gene expression by RT-qPCR, we 
used RNA isolated as described above, with three inde- 
pendent biological replicates, infected or not with X. 
fastidiosa. These RNAs were used for the cDNA synthe- 
sis according to the instructions of the Thermo Scientific 
for the RevertAid H Minus First Strand cDNA Synthesis 
Kit. After synthesis, the cDNAs were diluted at 1:25 and 
used in RT-qPCR. The evaluations were performed on 
an ABI Prism 7500 Sequence Detector System (Applied 
Biosystems, Foster City, CA, USA) using absolute quan- 
tification analysis. The detection of PCR products was 
measured by monitoring the increase in fluorescence 
emitted by SYBR green marker. For all amplifications 
performed in RT-qPCR, we produced dissociation curves 
to check for nonspecific amplification resulting from 
possible contamination. 

The analysis and normalization of gene expression 
were performed in the Genex software (version 5.0.1.5; 
www.multid.se), using the efficiencies and Cqs generated 
in the Miner software, and transformed into non- 
standardized data (Q). Uninfected samples were used as 
calibrators for each genotype evaluated. We used two 
endogenous genes for the normalization of the data. 

Additional files 



Competing interests 

The authors declare that they have no competing interests. 
Authors' contributions 

MAM and AAS planned and supervised the study. CMR and AAS contributed to 
the design and execution of the experiments. CMR conducted inoculation, 
collection, processing of samples, functional analysis of differentially expressed 
genes and drafted the manuscript. LTKand MAT contributed to the RNA-seq 
analysis, categorization and annotations of differentially expressed transcripts. 
CMR, AAS, MAT and MAM contributed to the interpretation of the data and 
provided intellectual input. MAT, AAS, LTK and MAM revised the manuscript. 
All authors read and approved the final manuscript. 

Acknowledgements 

The present study was conducted with the financial support of Fundagao de 
Amparo a Pesquisa do Estado de Sao Paulo (FAPESP) (Process n. 2007/07944-3), 
which provided a scholarship to CMR and Institute Nacional de Ciencia e 
Tecnologia (INCT) de Genomica para Melhoramento de Citros (Process n. 
573848/2008-4), with the financial support of FAPESP and Conselho Nacional 
de Pesquisa de Desenvolvimento (CNPq). MAM, AAS, and MAT are CNPq 
fellows. 

Author details 

^epartamento de Biotecnologia, Centro APTA Citros Sylvio Moreira, CP4, 
Cordeiropolis, SP 13490-970, Brazil. 2 Departamento de Genetica e Evolugao, 
Universidade Federal de Sao Carlos UFSCAR, CP676, Sao Carlos, SP 
13565-905, Brazil. 



Additional file 1: Detection of X. fastidiosa in plants of Ponkan 
mandarin and Pera sweet orange by RT-qPCR. DNA samples were 
prepared from xylem tissue from Ponkan mandarin and Pera sweet 
orange, after infection with X. fastidiosa or not (control) for one day, with 
five biological replicates for each species and their respective controls. 
(1 -5c) Pera sweet orange control; (6-1 0c) Ponkan mandarin control; (1-5d) 
Pera sweet orange with X. fastidiosa; (6-1 Od) Ponkan mandarin infected 
with X. fastidiosa; (c +) Positive control of a plant with CVC. Y-axis 
represents cycle quantification (Cq) as determined by RT-qPCR. 

Additional file 2: Differentially expressed genes in Ponkan 
mandarin after infection byX. fastidiosa (P< 0.001). Up and 

downregulated genes in Ponkan mandarin 24 h after infection with 
X. fastidiosa compared with the control. 

Additional file 3: Differentially expressed genes involved in the cell 
wall and secondary metabolism, according to MapMan results. Up 

and downregulated genes in Ponkan mandarin 24 h after infection with 
X. fastidiosa compared with the control that are related to secondary 
metabolism and the cell wall. 

Additional file 4: Differentially expressed genes involved in biotic 
stress, according to MapMan results. Up and downregulated genes 
related to biotic stress in tangerine Ponkan mandarin 24 h after infection 
with X. fastidiosa. 

Additional file 5: Differentially expressed genes involved in 
systemic ubiquitin dependent degradation, according to MapMan 
results. Up and downregulated genes related to ubiquitin dependent 
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